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O 1 ABSTRACT 

GO Aims. We present the results of a search for pulsed very-high-energy (VHE) y-ray emission from young pulsars using data taken with the 

H.E.S.S. imaging atmospheric Cherenkov telescope system. 
•> Methods. Data on eleven pulsars, selected according to their spin-down luminosity relative to distance, are searched for y-ray signals with 
•rH . periodicity at the respective pulsar spin period. Special analysis efforts were made to improve the sensitivity in the 100 GeV y-ray energy 
' domain in an attempt to reduce the gap between satellite and ground-based y-ray instruments. 

, Results. No significant evidence for pulsed emission is found in any data set. Differential upper limits on pulsed energy flux are determined for 
. . . i all selected pulsars in the approximate y-ray energy range between 100 GeV and 50 TeV, using different limit determination methods, testing a 
wide range of possible pulsar light curves and energy spectra. 

Conclusions. The limits derived here imply that the magnetospheric VHE y-ray production efficiency in young pulsars is less than 10~ 4 of 
the pulsar spin-down luminosity, requiring spectral turnovers for the high-energy emission of four established y-ray pulsars, and constrain the 
inverse Compton radiation component predicted by several outer gap models. 

Key words. Gamma-rays: observations - Pulsars: individual: Crab pulsar, Vela pulsar, PSRB1259-63, PSRB1509-58, PSRB1706-44 
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1. Introduction 

Rotating neutron stars are known to convert a significant part 
of their rotational energy into radiation originating from within 
the magnetosphere. This emission is observable as a periodic 
signal at the neutron star rotation frequency (the pulsar phe- 
nomenon). For many of the known young and energetic pul- 
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sars, the emitted luminosity p eaks at X-ray or y-ray ener 



gies (see iThompson et al.lll999l and references therein), usu 



ally attributed to curvature radiation of accelerated electrons 
in the strong magnetic fields pervading the pulsar magneto- 
sphere. The luminosity of the pulsed high-energy emission was 
found to correlate significantly with t he energy loss rate of the 
pulsar, i.e. its spi n-down power E (IDermer & Sturnerl 1 1 994 ; 



Cheng et al.l l2004). which relates to the magnetic field at the 
neutron star surface, fi sur f , approximately as E oc with 
Q. as the angular frequency of the pulsar rotatio n. For most 



of th e pulsars with esta blished y-ray emission ( Fie rro et al 
1 1 9981: iKuiper et al. 1999b . there is evidence for a turnover in 



the pulsed spectrum at a critical energy E c in the sub-GeV to 
lOGeV range. Note that no such turno ver has been found up t o 
lOGeV in the case of PSRB 1706-44 dThompson et alJll996h . 

Currently, two different scenarios for magnetospheric y- 
ray emission are discussed, placing the emission regions 
either near the m agnetic poles o f the neutron star {po- 
lar cap model, Sturrockl 197 it Ru derm an & Sutherland 



1971 iDaughertv & Harding! 1 1 996t iRudak & Dvksl 1 1 999t 



Dyks & Rudak 2004) or near the null surface in the outer 



magne tos phere of the pulsar (outer sap m o del. iCh eng et al. 
1986aH IChTang&Romanil Il992l ll994t iRomanil 1 199a 



Zhang & Chen g 1997; Cheng et al . 2000). Both models predict 



a cutoff in the curvature radiation spectrum at y-ray energies of 
the order of GeV up to several tens of GeV. However, the cut- 
off is of a different physical nature in the two models: while in 
the case of the polar cap model, the maximum energy of cur- 
vature radiation photons escaping the magnetosphere is limited 
by pair production in the polar magnetic field, the maximum 
energy of the accelerated electrons determines the maximum 
photon energy in the outer gap scenario. This results in differ- 
ent predicted light-curves, spectral shapes, and cutoff energies 
E c of the emitted high-energy radiation in the two models. Note 
that the predictions of outer gap models generally depend more 
strongly on model parameters, such as the inclination between 
the magnetic and rotational axes and the viewing angle to the 
observer, than is the case for polar cap models. Additionally, 
in some outer gap model calculations, a spectral component 
in the TeV range due to inverse Compton (IC) up-scattering 
of soft ambie nt seed photons by the accelerated electrons is 
predicted (e.g. Hirotani & Shibata 2001b; Hirotani et al. 2003; 



(e.g . Hir 

Takata et al.l2006l) . The boosted TeV photons, the flux of which 
is determined by the seed photon density, are more likely to es- 
cape due to the lower magnetic field strength in the outer mag- 
netosphere compared to the polar cap regions. 

Bridging the space between the distinct polar cap an d outer 
gap m odels, particle acceleration within the slot gap (lAronsI 



1983), a narrower region between the last open magnetic field 



line and the pair-plasma region of the polar cap, has been dis- 
cussed in the context of high-energy radiation. However, so far 
slot gap calculations have been used p rimarily to model de 



tails of pulsed high-energy l ight-curves (Dvks& Rudak 2003 



iMuslimov & Harding! 120031) without reproducing the corre- 
sponding energy spectra in the GeV range or predicting spectra 
for TeV energies. 

Pulsars with established y-ray emission components have 
been subject to intensive searches for pulsed very-high-energy 



(VHE, energies above ~ 100 GeV) y-ray emission by ground 
based instruments. Up to now, no evide nce for pulsed emission 
has been found in the s e observations dYoshikoshi et"aL 1997 



Chadwick et all l2000l; 



2002; Lessard et al 



Ide Naurois et aL _ _. _ - 
2000t lAharonian et alJl2004ah . and upper limits on the pulsed 
VHE y-ray flux have been derived under various assumptions 
on the characteristics of the pulsed emission. However, the IC 
component predicted by outer gap models has not yet been sig- 
nificantly constrained. 

The superior sensitivity of the H.E .S.S. detec tor with re- 
spect to previous instruments (see e.g. lOng 2006), combined 
with its Southern Hemisphere location, puts the predicted 
pulsed IC component from outer gap models within reach of 
testability, motivating a new search for pulsed VHE emis- 
sion from the known y-ray pulsars. A significant amount of 
H.E.S.S. data were accum ulated in a survey of the inner Galaxy 
(lAharonian et al. 2006d), resulting in the detection of numer- 
ous new sources of VHE y-rays, many of which are apparently 
associated with pulsars. Thus, in addition to the known y-ray 
pulsars, these data can be used to search for pulsed emission 
from promising young, energetic pulsars. 

The paper is structured as follows: in Section[2]the selection 
criteria of the studied pulsar sample are explained, Section|3]in- 
troduces the H.E.S.S. instrument and data analysis, in Section|4] 
the methods employed in the search for variability and the flux 
limit determination are explained, and finally, the results are 
presented and discussed in Sections |5]and|6] respectively. 

2. Pulsar Candidate Selection 

The prime candidates for the search for VHE y-ray emission 
are the pulsars with established y-ray emission at energies be- 
low ~ 10 GeV which have been detected by CGRO instruments. 
Four of them, PSRB0531+21 (Crab pulsar), PSRB0833-45 
(Vela pulsar), PSRB1706-44, and PSRB1509-58 were the 
subject of pointed H.E.S.S. observations and many hours of ob- 
servations were obtained for each of them in the years 2002 to 
2005. In addition, a significant number of known pulsars were 
observed as part of the H.E.S.S. survey of the Galactic Plane in 
2004 and 2005, or lie within the field of view of observations 
of other H.E.S.S. targets. The location of the H.E.S.S. site in 
the southern hemisphere limits the declination of observations 
to < 30°. 

Apart from the known y-ray pulsars, other candidates for 
which H.E.S.S. data were available were selected from the 
ATNF pulsar catalogue (see Manchester et al. I I2005L and ref- 
erences therein) if their spin-down flux E/D 2 was greater than 
10 35 erg s _1 kpc 2 . Table Q] lists all candidates chosen along 
with selected measured and derived characteristics collected 
from the literature. 

Some peculiarities of individual pulsars are pointed out in 
the following. PSR J1747-2958, situated in a bow-shock pul- 
sar wind nebula (PWN ) is detected as a compact X-ray source 
dGaensler et al.l I2004T) . The only binary pulsar in the present 
sample, PSRB1259-63 orbits a Be-star companion. Its mag- 
netic field is interestingly low (~ 3x 10 11 G), suggesting a high 
cutoff energy in polar cap m odels. It has been disco vered as an 
emitter of transient X-rays (ICominsky et al.l 1 19941) and VHE 
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Table 1. The characteristics of the selected pulsars taken from Manchester et alj (J2005). Period, P, distance, D, spin-down age, 
spin-down luminosity, E, and the corresponding value for E/D 2 , and calculated magnetic field strength at the neutron star surface, 
B sur f, and the light cylinder, B L c, are listed. The last column shows the rank in E/D 2 within the ATNF catalogue. 



Pulsar name 


P 


D 


Age 


lo gl0 


(ergs-) l0g - 


' E/etg s 1 \ 


B sm f 




Rank 




PSR 


[ms] 


[kpc] 


[kyears] 




\ £> 2 /k P c 2 j 


[10 n G] 


[10 4 G] 


E/D 2 


B0531+21 


J0534+2200 


33.1 


2 


1.24 


38.7 


38.1 


37.8 


98.0 


1 


B0833-45 


J0835-4510 


89.3 


0.29 


11.3 


36.8 


37.9 


33.8 


4.45 


2 


B 1706-44 


J1709-4429 


102 


1.8 


17.5 


36.5 


36.0 


31.2 


2.72 


6 


B 1509-58 


J1513-5908 


151 


4.4 


1.55 


37.3 


36.0 


154 


4.22 


7 




J1747-2958 


98.8 


2.5 


25.5 


36.4 


35.6 


24.9 


2.42 


13 


B 1259-63 


J1302-6350 


47.8 


1.5 


332 


35.9 


35.5 


3.3 


2.87 


15 




J1811-1925 


64.7 


5 


23.3 


36.8 


35.4 


17.1 


5.92 


18 




J1524-5625 


78.2 


3.8 


31.8 


36.5 


35.3 


17.7 


3.46 


19 




J1420-6048 


68.2 


7.7 


13 


37.0 


35.3 


24.1 


7.13 


22 




J1826-1334 


101 


4.1 


21.4 


36.4 


35.2 


27.9 


2.51 


23 




J1801-2451 


125 


4.6 


15.5 


36.4 


35.1 


40.4 


1.95 


30 



y-rays dAharonian et al. 2005 ). however, pulsed high-energy 
emission has not so far been detected. PSR J 1420-6048 is em- 
bedded in a PW N and identified with a compact X-ray source 
dNg et al.l 120051) within the Kookaburra com plex, from which 
VHE y-ray emission was recently detected dAharonian et al 



2006a). PSR J181 1-1925, for whi ch pulsed emission has been 
detected up to hard X-rays dGavriil et alJ2004h . is located in the 
composite supernova remnant Gl 1.2-0.3. All of the selected 
pulsars are visible at radio wavelengths with the exception of 
PSR J181 1-1925 which is only seen in X-rays. 

All available H.E.S.S. data were searched for the respective 
pulsar to be in the field of view with a maximum angular dis- 
tance of 6 — 2° between the pulsar position and the H.E.S.S. 
pointing direction. 

3. VHE Gamma-Ray Detection and Data Analysis 

3.1. The H.E.S.S. Detector 



The High Energy Stereoscopic System (H.E.S.S., Hof mann 
2003I) . an array of imaging atmospheric Cherenkov telescopes 
located in Namibia at 23°16'18"S 16°30'00"E and 1800m 
a.s.l., detects cosmic VHE y-rays by imaging the Cherenkov 
emission of their air showers in the atmosphere using opti- 
cal telescopes. Each telescope has a tess ellated spherical mir 



ror with 13m diameter and 107 m area ( Bernloh r et al.ll2003 



Cornil s et al. 2003) and is equipped with a camera of 960 
0.16°-photomul tiplier tubes providi ng a total field of view of 
5° in diameter ( Vincent et alJ 2003 ). During the stereoscopic 
observations, an array trigger requires the simultaneous de- 
tectio n of air showers by several telescopes at the hardware 
level ( Funk et all 20041) . allowing a suppression of background 
events and the assignment of GPS timestamps to each triggered 
air shower, without the need of an offline synchronisation of 
the individual telescopes. After digitisation, t he event data are 



stored on disk by the data acquisition system ( Borgmeier et al 
20031) . 

Fo llowing the cal i bration of the recorded air shower 



data dAharonian et al.1 2004b), each telescope image was 



the shower geometry providing an angular resolution of ~ 0.1° 
for individual y-rays. The y-ray energy was estimated from the 
image intensity and the shower geometry with a typical resolu- 
tion of ~ 15%. In order to reject the vast background of cosmic- 
ray showers, y-ray candidates are selected using cuts on im- 
age shape scaled with their expectation values obtained from 
Monte Carl o simulations. De t ails of t he analysis technique can 



parametrised by its centre of gravity and second moments 
dffillasl 985) followed by the stereoscopic reconstruction of 



be found in lAharonian et al.l (l2006cl) . The selection cuts used 
for this analysis were optimised on simulations of a y-ray point 
source with 10% of the flux of, and the same spectral shape as, 
the VHE y-ray emission from the Crab Nebula. These cuts, in 
the following called standard, allow the detection of sources 
with 1 % Crab Nebula flux at energies above ~ 200 Ge V with 
a significance of 5 <x within 25 hours at a zenith angle of 20°. 
Note that the threshold energy above which y-ray sources can 
be detected increases with the zenith angle of observation. 



3.2. Low Energy Optimisation 

Since observational data indicate steep cut-offs in high-energy 
(GeV) y-rays, special low-energy cuts have been applied to re- 
duce the gap in observational coverage between satellite and 
ground based y-ray observations of young pulsars. Table|2]lists 
the value s of the two samples of se lection cuts used in this anal- 
ysis. See lAharonian et al. (2006c) for a detailed description o: ' 
the individual cuts. The low energy cuts were optimised on sim- 
ulated y-rays to yield the best sensitivity for steep spectrum 
sources (dN/dE oc E _r ,r = 5) with energies below 500GeV. 
Figure Q] shows the rate of simulated y-ray events for a Crab- 
like source spectrum (r = 2.6) passing the standard (dark grey) 
and low energy (light grey) selection cuts as a function of y-ray 
energy for a typical zenith angle of 20°. The peak of the distri- 
bution dN/dE is commonly defined as the threshold energy. 
The low energy cuts result in an energy threshold lowered by 
a factor of 1.5 with respect to the standard cuts. Although the 
number of selected events is higher by 50%, the background re- 
jection efficiency of the low energy cuts is worse by more than 
a factor of 3. Thus, the extension of the detectable y-ray energy 
range is achieved only at the expense of decreased sensitivity, 
especially at higher energies. 
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10 _ 10 

Energy [TeV] 

Fig. 1. Event rate of simulated y-rays as a function of y-ray 
energy for two different sets of selection cuts. The dark grey 
histogram corresponds to the standard and the light grey his- 
togram to the low energy cuts. The vertical arrows roughly rep- 
resent the corresponding threshold energy. 

Table 2. Selection cuts used for the analysis of air shower 
events. The image amplitude (Size) is given in photo-electrons 
(p.e.), E lmx corresponds to the maximum reconstructed event 
energy, MRSW and MRSL represent the mean scaled reduced 
width and length parameters, and #;? ut is the maximum squared 
angular distance with respect to the target position. The tail cuts 
represent the image clea ning thresholds applied. For details see 
Aharonian et al . (2006c). 



Cuts 


Size 


F 

^max 


MRSW 


MRSL 


c 


Tail cuts 




[p.e.J 


TeV 






Ldeg 2 ] 


[p.e.] 


Standard 


>80 


100 


[-2,0.9] 


[-2,2] 


0.0125 


5/10 


Low energy 


>40 


0.5 


[-2, 0.75J 


[-2,1.4] 


0.024 


4/7 



In order to quantify possible systematic effects of the 
low energy cuts on the energy spectrum determination, data 
from observations of the Crab Nebula were used to recon- 
struct the differential energy spectrum for both sets of cuts. 
In Figure [2] the open points show the spectrum of unpulsed 
y-rays determined with the standard cuts. The solid line indi- 
cates the best-fit power l aw with exponential cutoff taken from 
Aharonian et al.1 (|2P06c), Table 5. The full points result from 
the analysis with low energy cuts, showing minor systematic 
differences in the absolute flux normalisation (<10%). 



3.3. Point Source Analysis 

Observation runs are selected according to standard quality se- 
lection criteria, ensuring sta ble detector and at mospheric con- 
ditions as described in lAharonian et aL ( 2006c ). The data were 
analysed using the standard method for the search for point 
sources without considering any time variability. The number 
of events in a circular region around the source position with ra- 
dius cut is compared to the number of counts in several back- 
ground regions of the same size, arranged on a circle around 




10' 

Energy [TeV] 

Fig. 2. Differential energy spectrum of unpulsed VHE y-rays 
from the Crab Nebula obtained with the different sets of selec- 
tion cuts. The open points correspond to the results of the stan- 
dard and the full points to the low energy cut analysis (but with 
£ mal = 100 TeV as for the standard cuts). The solid line repre- 
sents th e spectrum parametrisation taken from I Aharonian etaD 
(l2006ch . 



the centre of the field of view with radius o g. For each run, 
the background normalisation is determined by the number of 
background regions used. For data sets comprised of several 
different offsets, as is the case for most of the selected pulsars, 
the background normalisation for each run was weighted by 
its relative acceptance for a given offset. This acceptance was 
derived from data of H.E.S.S. observations without significant 
y-ray emission within the field of view (see also iBerge et all 
2006h . 



4. Periodicity Analysis 

In order to test for pulsed y-ray emission at the pulsar posi- 
tion, the timestamps of each recorded shower passing selection 
cuts were transformed from the observer's frame into the pul- 
sar frame and then folded with the pulsar spin period. The re- 
sulting unbinned distribution of pulsar phases corresponding to 
each shower event was tested for variability. In cases where no 
significant periodicity was found, flux limits were determined 
under various assumptions concerning pulse shape and posi- 
tion. 



4.1. Timing Corrections 

The time of arrival (TO A) for each shower event recorded by 
the tel escope array is d etermined with a central trigger GPS 
clock dFunk et alJbOoH) . achieving an accuracy of better than 
2 yus. Each TOA fxoA was transformed to the pulsar frame using 
the following corrections: 



t = f T0 A + Af dock + Af SSB + At; 



Shapiro 



+ At] 



Einstein 



At] 



binary ■ 



(l) 



Here, Af c i oc k represents the transformation from GPS time 
to International Atomic Time, AfssB the classical transfor- 
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mation to the solar system barycentre using the solar sys- 
tem ephemeride s DE200 from the Jet Propulsion Laboratory 
dStandishlll982l) and Af shapil0 and Af Einstein are general rela- 
tivistic corrections for t he gravitational influence of the sun 
(Bac ker & Hel lings 119861) and the redshi ft induced by all ma- 
jor bodies in the solar system dFairhead et alJll988l) . respec- 
tively. The correction term Aginary applies only for pulsars in 
binary systems, i.e. PSRB 1259-63, a nd corrects for the pulsar 
motio n in its binary orbit according to Bla ndford & Teukolsky 
dl976h . The correction software used for this analysis was 
cross-che cked with the com monly used pulsar timing utility 



TEMPO (Tay lor et alj 12000). Furthermore, the timing capa- 
bility of H.E.S.S. was demonstrated by measurin g the optical 



light-curve of the Crab pulsar with high precision ( H inton et al 
2006). 



4.2. Tests for Pulsed Emission 

The corrected TOA was subsequently folded with the timing 
parameters of the target pulsar (pulsar ephemeris) obtained 
from radio or X-ray observations. For each event TOA f, the 
pulsar phase <p was determined according to 



4>{t) = NAt) - LJV r (0J 



(2) 



where |_. ■ -J denotes the floor function, and with iV r as the num- 
ber of pulsar rotations 

1 , 1 , 

N T (t) = M(f0) + V\ t0 (t - t Q ) + -H„(f - to) 2 + -H„(f - toY (3) 



and v|, are the 



where the parameters fo, N T (to), v|, , 
reference time, the corresponding number of rotations, fre- 
quency, and first and second frequency derivative, respectively, 
which altogether represent the timing parameters of the pulsar 
ephemeris. Table [3] lists the references for the ephemeris used 
for each pulsar along with the respective reference times to de- 
pending on the number of different parameter sets. 



Table 3. Pulsar ephemerides. The first and last reference times 
(fo) of the timing solutions used in this analysis are given. 
For a complete list of the solutions a nd the correspon ding 



pulsar parameters see the r efer ences: 1 : iLvne et al 
Manchester, R. et al.1 d2005l) . 3: ICamilo et al.ld2002h 



(2005J), 2: 



Pulsar 


t [MJD] 


References 


B0531+21 


52929,..., 53446 


(1) 


B0833-45 


53 123,..., 53469 


(2) 


B 1259-63 


53108,..., 53386 


(2) 


J1420-6048 


53156, 53493 


(2) 


B1509-58 


53111, ...,53385 


(2) 


J1524-5625 


53650 


(2) 


B 1706-44 


52773, 52915 


(2) 


J1826-1334 


53130, 53449 


(2) 


J1747-2958 


52613* 


(2) 


J1801-2451 


52503* 


(3) 


* Valid timing 


information was not available for 


all data within the H.E.S.S. data set. 





For the X-ray pulsar PSRJ1811-1 925 and the radio pulsars 
PSR J1801-2451 and PSR J1747-2958 there were no contem- 
poraneous observations available from which the pulsar tim- 
ing parameters could have been deduced. For these pulsars, the 
most recently published timing solution was used, assuming 
that the pulsar phase determined with the known parameters 
is accurate enough to reveal the periodicity of a potential sig- 
nal. In particular, this assumes that no glitches or other major 
frequency changes occurred within the time range between the 
H.E.S.S. observations and the measurement of the timing pa- 
rameters. 

The unbinned phase distribution of events was tested for 
significant periodicity, i.e. for a significant deviation from a flat 
distribution, using several statistical tests. The Pearson ^- 2 -test 
tests for a constant distribution of phases in a number of phase 
bins, which was set to 20 in order to give acce ptable coverage 
for an y of the considered data sets. The Z 2 -test (I Buccheri et al 



1983b checks for a combination of fundamental sine and co- 
sine harmonics up to the ra-th order and is independent of any 
binning. Here, the Z\- and Z| -tests for single an d double peak 



profiles are employed. Additionally the //-test dde Jager et al 
1989), a more powerful generalisation of Z„,, was performed 
to increase the sensitivity to unknown pulse profiles. Finally, 
a rather different approach which checks for maximum de- 
viations from the uniform phase distribution, the Kuiper-test 
dJetsu & Peltll 19961) . was used. 

In order to test for spurious signals induced by possible sys- 
tematic effects in the H.E.S.S. event data, the tests were also 
applied to the phase series of background events, i.e. events 
from sky regions used in the background determination for the 
unpulsed y-ray excess (cf. Sec. 13.31 ). 

As an example, the distribution of event phases from ob- 
servations of the Vela Pulsar (PSRB0833-45) is shown in 
Fig. [3] obtained using the standard cuts. The difference be- 
tween on and off results from the known y-ray excess from 
HESS J0835-456 at the position of the pulsar (see Tabled). 
No significant deviation from uniformity was found with any 
of the statistical tests for pulsations. 

4.3. Flux Limit Determination 

Several methods were applied to obtain limits on the y-ray flux 
from the selected pulsars. They differ in the assumptions made 
concerning the characteristics of the pulsed emission, in par- 
ticular concerning the form of the light-curve and the phase 
position. 

In the most conservative approach, without any assump- 
tions on the temporal structure or spectrum of a possible signal, 
the energy flux E 2 dN/dE is calculated using the standard anal- 
ysis as described in Sec. 13.31 This method considers all events 
in the direction of the pulsar, i.e. in the on-region, ignoring any 
possible variations and is referred to as the unpulsed method. 
The flux is calculated event-wise for the on- and off- (back- 
ground) regions within the same field of view separately, accu- 
mulated in bins of energy, and sub sequently subtracted to yield 
the excess flux dAharonian et al.ll2005cl) . Note that due to the 
finite width of the energy bins, the reconstructed flux is weakly 



6 



F. Aharonian et al.: H.E.S.S. Search for Pulsed VHE y-rays from Young Pulsars 




On Region: 

P(H) = 0.92, P(Z?) = 0.9, P(Zl) 
P(Kuiper) = 0.88, P X 2= 0.79 
Pulsed Fraction (Zf ) < 0.058 (99% CL) 

Off Regions: 

P(H) = 0.84, P(Z?) = 0.81 , P(Zl) = 0.76 



P(Kuiper) = 0.68, Py2 
Pulsed Fraction (Z|) 



= 0.49 
: 0.046 (99% CL) 



H.E.S.S. 

live time = 1 5 h 

Vela Pulsar 




0.9 1 

Pulsar Phase 

Fig. 3. Upper plot: Distribution of event phases for the Vela 
pulsar (PSRB0833-45). The points represent the events in the 
on-region at the pulsar position and the histogram the nor- 
malised off-region events. The solid and dashed lines represent 
the x 1 fits to a constant for the on- and off-events, respectively. 
For both regions the probabilities for being consistent with a 
uniform distribution according to the statistical tests on pulsa- 
tions (see text) are listed along with the pulsed fraction limit as 
introduced in Sec. 14. 31 Lower plot: Phase distribution for y-rays 
with energies betw een 2 and lOGeV as measured by EGRET 
(iFierro et al. I ll998h . The dotted vertical lines denote the edges 
of the on-pulse regions. 



dependent on the y-ray energy spectrum for which the detector 
efficiency was calculated. Finally, the flux obtained in this way 
is convert ed into a 99% confidence le vel upper limit using the 
method of iFeldman & Cousins! U.998). As some of the pulsars 
in this sample show significant unpulsed y-ray emission, the 
flux limits derived can be considerably higher than the nominal 
H.E.S.S. y-ray flux sensitivity since the unpulsed y-ray emis- 
sion is an additional background to the pulsed emission. 

The second method takes advantage of the pulse position 
and shape as measured in other energy domains, assuming sim- 
ilar characteristics for the hypothetical VHE y-ray emission. 
In particular for the known y-ray pulsars this approach seems 
reasonable, since the y-ray emission mechanism is expected 
to be similar at VHE y-ray energies, especially for photons 
originating from curvature radiation. For each pulsar, a known 
pulse profile was chosen according to the energy band of the 
emission, taking the highest photon energies. For these pulse 
profiles, an on- and an off-phase domain were assigned (pn- 
off-pulse method), representing the pulsar phases for high and 
low (or no) emission, respectively. Table [4] lists the used on- 
phase intervals and their corresponding energy bands for all 
pulsars. The phase intervals are given relative to the reference 
phase which is coincident with the highest peak of the respec- 



tive radio profile. For example for the Vela pulsar, the on- and 
off-phase domains were chosen based o n the light-curve m ea- 
sured using EGRET at GeV energies dFierro et al. 1998 ). as 
indicated by the dotted lines in Fig. [3] Note that any possi- 
ble contribution from inter-pulse domains, the phase range be- 
tween the on-pulse domains, to the pulsed signal are neglected 
in this method. For the pulsars without established y-ray emis- 
sion, the selected phase region always covers 20% of the pro- 
file around the reference phase. The VHE y-ray flux was then 
calculated as for the unpulsed method, but using the off-phase 
events in the direction of the pulsar for background estima- 
tion, instead of the events in the off-regions within the field of 
view. If there was no significant unpulsed signal visible within 
a data set and no known VHE y-ray source coincident with 
thetargetpulsar,i.e.forPSRJ1524-5625,PSRB1706-44, and 
J 180 1-2451, the off-event statistics were increased by includ- 
ing those events which were in the off-regions of the field of 
view and additionally had phases falling in the assigned off- 
phase domain. In these cases the background normalisation was 
adjusted accordingly. The on-off-pulse method is less reason- 



Table 4. Pulse phase regions used in the on-off-pulse method 
Given are the phase intervals considered as on phase to- 
gether with the energ y band and the c orre sponding ref 
erence. References: 1: IFierro et al. (1998), 2: iKramer et al 
d2003h. 3:lKuiperetalJ dl999h. 4: iThompson et all dl996l). 5 
ICamilo et alJ d2002l) . 7: lTorii et al.l (11997^ " 



Pulsar 


On pulse phases 


Energy 


Reference 


B0531+21 


[-0.06, 0.04], [0.32, 0.43] 


y-ray 


(1) 


B0833-45 


[0.1,0.14],[0.5,0.57] 


y-ray 


(1) 


B1259-63 


[-0.1,0.1] 


Radio 


(2) 


J1420-6048 


[-0.1,0.1] 


Radio 


(2) 


B 1509-58 


[0.15,0.65] 


y-ray 


(3) 


J1524-5625 


[-0.1,0.1] 


Radio 


(2) 


B 1706-44 


[0.24,0.5] 


y-ray 


(4) 


J1826-1334 


[-0.1,0.1] 


Radio 


(2) 



able for those pulsars for which pulsed emission is known only 
at low energies, especially at radio wavelengths, since the y-ray 
production mechanism may be completely unrelated to that of 
the low energy emission. 

Therefore, a third approach was also applied, that does not 
rely on knowledge of the phase position of the pulse, but only 
assumes a certain pulse shape. Here, a combination of first and 
second order harmonics was chosen, similar to the Z|-test. A 
99% confidence level upper limit on the fraction, p<Si 1, of 
pulsed emission present in the series of event phases for events 
with shower directions from within the on-region around the 
pulsar position was c alculated for a given energy bin accord- 



ing to Brazierf ( 1994 ). The resulting event fraction was con- 



verted into a flux using the effective collection area of the de- 
tector and assuming a spectral distribution in energy following 
a power law with photon index F = 2.7, roughly matching the 
shape of the cosmic ray background s pectrum. More details o n 
the flux calculations can be found in dAharonian et al.l l2006c). 
This method, in the following referred to as the pulsed frac- 
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tion method, is less constraining than the on-off -pulse method, 
since no background subtraction can be performed and thus the 
post-cut cosmic ray background is treated as an unpulsed com- 
ponent in addition to any unpulsed VHE y-ray excess. 

Altogether, the pulsed fraction method was considered to 
be the most robust approach since the structure of the VHE 
y-ray light curve remains unknown a priori. Nevertheless, the 
on-off-pulse method was still applied to ensure comparability 
with results from prior searches for pulsed VHE y-ray emission 
where very similar techniques have been used. For the three 
pulsars for which valid timing parameters were not available 
for the H.E.S.S. data, PSRJ1747-2958, PSRJ1801-2451 and 
PSRJ181 1-1925, flux limits have been derived only using the 
unpulsed method. 



5. Results 

5.1. Point Source Analysis 

Table[5]shows the main characteristics of the data sets for each 
pulsar along with the results of the point source analysis for un- 
pulsed emission. Additionally, any H.E.S.S. VHE y-ray source 
close to the corresponding pulsar position is listed whenever a 
significant y-ray excess was detected which can be attributed to 
this known source. Note that the excess given may differ signif- 
icantly from that given in the respective reference, since in all 
cases, with the exception of the Crab nebula, the VHE emission 
region is both extended and offset from the pulsar position. 

Three of the four pulsars not coincident with known 
VHE y-ray sources - PSR J1524-5625, PSRB1706-44, 
and PSR J 180 1-2451 - show no significant signal, but for 
PSR J181 1-1925 an excess with 3.8 cr significance is seen. 
After taking into account all four trials for these target posi- 
tions, the significance is reduced to 3.5 cr. 

Since the H.E.S.S. energy threshold increases with the 
zenith angle of observation, the low energy selection cuts in- 
troduced in Sec. 13.21 effectively reject all data taken at zenith 
angles higher than 50° since the effective collection area of the 
instrument vanishes below 500 GeV y-ray energy for these el- 
evations. Therefore the effective exposure time for the low en- 
ergy analysis is reduced for some data sets compared to that of 
the standard cuts. 

For PSR B 1706-44, only archival data from the early com- 
missioning phase of H.E.S.S. from observations with two 
telescopes and without an array level trigger are available. 
Therefore, event timing was determined from the GPS times- 
tamps of the individual cameras. Additionally, the low energy 
analysis was not performed. 

In the special case of PS RB1259-63, the data set used 
here differs from that used in I Aharonian et alj (12005I) . Newer 
data from observations in 2005 are included and data taken 
before MJD 53135 are excluded. Before this date reliable ra- 
dio ephemerides could not be provided as the radio pulsar was 
eclipsed by its binary companion. Additionally, the unpulsed 
y-ray emission, presumably arising from the interaction of the 
pulsar wind and the stellar outflow, would result in a strong 
background for any pulsed emission component. 
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On Region: 

P(H) = 0.0031, P(Zf) = 0.072, P(Z|) = 0.007 
P(Kuiper) = 0.0049, P Z 2= 0.0093 
Pulsed Fraction (Zl)<0.21 

Off Regions: low energy cuts 

P(H) = 0.64, P(Z?) = 0.57, P(Z|) = 0.67 live time = 1 2 h 

PSR J 1420-6048 



P(Kuiper) = 0.39, Py2= 0.082 
Pulsed Fraction (Zf) < 0.034 
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Fig. 4. Distribution of event phases for PSR J 1420-6048 using 
the low energy cut analysis for on- and off-regions (see also 
Fig. |3). Note that there is significant unpulsed emission at the 
pulsar position resulting in a difference between the mean num- 
ber of events for the on and off-regions shown by the solid and 
dashed lines respectively. 



5.2. Search for Pulsations 

Table [6] shows the results of the statistical tests applied to 
the phase distributions of events in the on-region around each 
pulsar position. The probabilities for the hypothesis of a flat 
phase distribution, i.e. the probability that no pulsed signal 
exists, are given for the two sets of selection cuts applied. It 
can be seen that these probabilities lie well above the one per 
cent level with a single exception. For the low energy analysis 
of PSR J1420-6048, the test probabilities seem significantly 
low, especially for the more powerful H- and Kuiper-tests. No 
such effect is seen for the respective probability of the back- 
ground distributions of any data set, which are all above 1%. In 
Figure |U the corresponding phase distribution is shown, again 
for both on and off-events. The off-event distribution is much 
more compatible with a flat distribution, with probabilities well 
above 10%. The probability PhIpsrji420-6048 = 3.1 X 10~ 3 of 
the H-test for the on-events corresponds to 3.0 cr. However, 
this probability has to be corrected for the number of trials 
for this test. Considering 19 trials for all data sets analysed 
with the low-energy (9 trials) and standard (10 trials) selec- 
tion cuts respectively, the corresponding significance for the 
applied tests lies between 2.0 cr (H-test) and 0.4 cr (Z^-test) 
with an average of 1 .5 cr. Note that the test results for the 
pulsars PSRJ1747-2958 and PSR J 180 1-2451 for which no 
valid timing information was available do not allow firm con- 
clusions about the presence of periodicity within the H.E.S.S. 
data with respect to the intrinsic pulsar rotation period. For 
PSR J181 1-1925, the tests were not performed since there was 
no valid timing information available for any period within sev- 
eral years of the H.E.S.S. observation time. 

As a cross-check for systematic effects possibly biasing the 
statistical tests, the ensemble distribution of the test statistic 
was checked against the hypothesis that no pulsed signal exists 
in any of the data sets. Figure shows the distribution of the 
test statistic H of the H-test for the on- and off-regions and for 
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Table 5. Results of the analysis of unpulsed y-ray emission for each pulsar data set using the standard cuts and assuming a point 
source at the pulsar position. Dead-time corrected observation time, f b S , mean offset of the pulsar direction with respect to the 
array pointing direction, (6 g), approximate mean energy threshold, (^ttesh)) number of excess counts, N excess , and corresponding 
significance, S, are given. In addition to the mean threshold for the standard cuts, the corresponding threshold for the low 
energy cuts is given in parentheses if applicable. If the pulsar is in the vicinity (angular distance < 1°) of, or associated with, 
a known VHE y-ray source, the source name is given along with the corresponding reference. Note that the given numbers 
can significantly differ fr om the published results for the known sources whic h ar e located at different po sitio ns and might be 



extende d. References: A: Aharonia n et al.l d2006cl). B: Aharonia n et al] (2006b), C: Aharonian et al. (2005), D: Aharonian et al 



d2006al) . E: lAharonian et al.l d2005bl) . F: lAharonian et al.l d2005dl) . G: lAharonian et alJ (l2006dl) . H: lAharonian et alJ J2005al) . 



Pulsar 


fobs 


(Ooff) 


(-^"(hresh) 


N 

1 ' excess 


5 


H.E.S.S. 


Reference 


PSR 


M 


[°] 


[GeV] 






catalogue source 




B0531+21 


25.4 


0.6 


510(310) 


9673 ± 59 


164.8 


Crab Nebula 


(A) 


B0833-45 


16.3 


0.6 


250 (170) 


97 ±40 


2.81 


HESSJ0835-456 


(B) 


B 1259-63 


55.8 


0.6 


380 (280) 


1081 ±65 


16.7 


HESSJ1303-638 


(C) 


J1420-6048 


13.0 


0.8 


350 (240) 


338 ±28 


12.07 


HESSJ1420-607 


(D) 


B 1509-58 


35.8 


0.6 


310(210) 


1853 ±58 


32.22 


HESSJ1514-591 


(E) 


J1524-5625 


15.5 


1.0 


330 (220) 


-8 ±28 


-0.30 






B 1706-44 


16.6* 


0.5 


255 


1 ±24 


0.3 




(F) 


J1747-2958 


103.3 


1.1 


200 (140) 


291±82 


3.62 


HESS J 1745-303 


(G) 


J1801-2451 


9.8 


1.2 


180 (130) 


-19 ±28 


-0.71 






J1811-1925 


14.2 


1.1 


210(150) 


116 ± 31 


3.8** 






J1826-1334 


62.8 


1.1 


210(140) 


683 ± 63 


11.2 


HESSJ1825-137 


(H) 



*Comprises 14.4 h of 2-telescope data and 2.2 h of data from the full 4-telescope array. 



See text. 



Table 6. Results of statistical tests for pulsed emission. Number of events, N on ,X 2 probability for a fit of the phase distribution 
to a constant, and probabilities for the event series being compatible with a flat phase distribution for the H-, Z^-, Z?-, and 
Kuiper-test (Pk) are shown. All quantities are given for both the low energy and standard cut analyses (see text). 



Pulsar 






Low energy cuts 










Standard cuts 






PSR 


A'on 


p * 


Ph 


z i 


P A 


Pk 


A'on 


p x- 


Ph 


P z i 


p zk 


Pk 


B0531+21 


8095 


0.99 


0.84 


0.81 


0.94 


0.97 


10622 


0.51 


0.56 


0.63 


0.25 


0.66 


B0833-45 


7480 


0.52 


0.87 


0.85 


0.37 


0.82 


1156 


0.79 


0.92 


0.90 


0.98 


0.88 


B 1259-63 


16176 


0.78 


0.92 


0.90 


0.65 


0.71 


4535 


0.29 


0.25 


0.18 


0.46 


0.23 


J1420-6048 


2228 


0.0093 


0.0031 


0.072 


0.007 


0.0049 


968 


0.67 


0.60 


0.53 


0.53 


0.62 


B 1509-58 


12481 


0.37 


0.81 


0.77 


0.70 


0.89 


4308 


0.048 


0.055 


0.027 


0.11 


0.04 


J 1524-5625 


2498 


0.78 


0.43 


0.35 


0.43 


0.39 


745 


0.87 


0.79 


0.75 


0.96 


0.97 


B 1706-44* 














391 


0.02 


0.82 


0.78 


0.85 




J1826-1334 


14497 


0.71 


0.46 


0.38 


0.57 


0.62 


4016 


0.46 


0.42 


0.34 


0.34 


0.42 


J1747-2958** 


23482 


0.65 


0.98 


0.97 


0.82 


0.95 


6340 


0.62 


0.62 


0.92 


0.99 


0.96 


J1801-2451** 


3230 


0.035 


0.22 


0.15 


0.42 


0.21 


723 


0.50 


0.15 


0.094 


0.24 


0.32 



*Subject to analysis limitations (archival data, see text). 



Valid pulsar timing information was not available for all data within the H.E.S.S. data set. 



the different sets of cuts. Also shown is the expected distribu- 
tion Nu(H) = No exp(-AH) for the case when no pulsed signal 
is present in any data set, with No depending on the number of 
tests and A - 0.4. In order to quantitatively evaluate the match 
between data and expectation, each of the distributions for on- 
and off-events and for the two sets of cuts was fitted by the 
respective expected distribution with a log-likelihood method, 
leaving No and A as free parameters, and the probability for be- 
ing compatible with the expected parameters was calculated. 
No significant deviation is found for any of the four distribu- 
tions with all probabilities above 20%. 



5.3. Flux Upper Limits for Pulsed Emission 

Differential energy flux upper limits (i.e. in E 2 dN/dE) were 
derived for each data set, for the energy range from thresh- 
old energies up to the highest energies for which significant 
event statistics were accumulated. In Fi gures loT-fTTlthe differen- 
tial flux limits calculated are shown separately for each pulsar. 
Different points are shown for the different limit determination 
methods (see Sec. 14. 3I >. The points above energies of 0.5 TeV 
correspond to the results for the standard selection cuts and 
the points below 0.5 TeV to these for the low energy cuts. The 
horizontal position of the points denotes the mean energy of 
the events in each energy bin while the bin edges are indicated 
by the horizontal error bars. In the case of the pulsed fraction 
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Fig. 5. Ensemble distributions of the H-test statistic for the se- 
lected pulsars of Table [6] and their corresponding background 
control samples. The results for the pulsar direction are shown 
as open and closed circles for the low energy and standard cut 
analysis, respectively. The distributions for the off-regions are 
displayed as grey filled and outlined histograms, respectively. 
The solid curve shows the expected distribution if no pulsed 
signals are present. 

method, the mean energy was calculated for the assumed power 
law energy spectrum. 

Note that the overall y-ray energy scale for the standard 
cuts has a systematic error estima ted to be 20% for a p ower law 
energy spectrum with F ~ 2.6 (I Aharonian et al. 2006cl) . domi- 
nated by the absolute shower energy calibration of the instru- 
ment. The impact of the systematic uncertainty in the energy re- 
construction method on the flux upper limit for the low-energy 
cuts was determined to be 30% for energies around the energy 
threshold of the respective data set, mainly affecting the given 
flux limit point for the lowest energy bin. 



6. Discussion 

For the four y-ray pulsars, these new results allow some general 
conclusions on high-energy emission mechanisms of pulsars to 
be drawn. The results are put in the context of model predic- 
tions from selected polar cap and outer gap scenarios. In the 
particular case of the outer gap scenario, only those model cal- 
culations were considered which make explicit predictions for 
the IC component at TeV energies. Where appropriate, simu- 
lated spectr a are shown which are calc ulated using simulatio n 
codes from iRudak & Dykl (1 19991) and iHirotani et all d2003h . 
Finally, the results can be compared with findings of other ex- 
periments in the TeV energy domain. 

Crab pulsar. The Crab pulsar has been studied by all major 
ground-based y-ray ins truments at e nergies between 60 GeV 
dde Naurois et"ai]l2002h and 50 TeV (lAharonian et al.ll2004ah . 
Figure|6]shows the new H.E.S.S. results along with y-ray mea- 
surements from EGRET and selected results in the VHE do- 
main. The southern location of H.E.S.S. allows only observa- 
tions at rather high zenith angles, prohibiting a deep exposure 
especially at low energy thresholds. Therefore the limits ob- 
tained using the Northern Hemisphere instruments HEGRA, 
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Fig. 6. H.E.S.S. energy flux limits (99% c.l.) for pulsed emis- 
sion of the Crab pulsar. The full circles and full squares corre- 
spond to the on-off-pulse and pulsed fraction limit determina- 
tion methods, respectively. Below energies of 0.5 TeV the re- 
sults were obtained with the low energy selection cuts, oth- 
erwise the standard cuts were used. Note that the threshold 
energy of the H.E.S.S. data lies at 3 10 GeV due to the ob- 
servation at high zenith angles imposed by the limited ac- 
cessibility of the northern sky. The crosses, open squares, 
and open diamonds de note limits from CELESTE (99 % c.l., 
I de Naurois eTal] 12002b , Whipple ( 99.9% c.l., ILessard et al 



Aharonian et al 



2004a) 



120001) . and HEGRA (99.865% c.l., 
respectively. The o pen triangles sho w the EGRET phase- 



averaged spectrum (iFierro et al.l 1 19981). The indicated polar 



cap curve was generated according to Rudak & Dyksl ( 1999) 
and the outer gap model curve taken from Hirotani & Shibata 
d2001al) . 



with more than 200 hou rs of observations, an d Whipple, with a 
lower energy threshold (ILessard et alJl2000i) . could not be im- 
proved upon. However, the HEGRA limits are derived only 
for a rather narrow pulse region between -0.06 and 0.04 
of the main pulse, whereas the Whipple and H.E.S.S. limits 
include both the main and secondary pulse regions. Taking 
the H.E.S.S. results for the on-off-pulse method and assum- 
ing a power-law spectrum with exponential cutoff dN/dE oc 
E~ r exp(-E I E c ), with T = 2.05, the cutoff energy is con- 
strained to be E c < 80 GeV, less constr aining than the flux limi t 
at -60 GeV reported by CELESTE dde Naurois et al.ll2002l) . 



Our new differential flux upper limits for the Crab pulsar 
constrain the IC compone nt in the TeV range predicted by 
dHirotani & Shibata 2001a), reinforcing the constraints placed 
by previous y-ray experiments, both at GeV and TeV energies. 
No constraints can be placed on the polar cap prediction, which 
exhibits a sharp cutoff at a few GeV. 

Vela pulsar. For Vela, the l imits set using the Durham 
Mark 6 telescope (IChadwick et al.l2o6oh and the CANGAROO 
(lYoshikoshi et al.lll997l) detector are improved by more than 
one order of magnitude in flux and reach down to ~ 100 GeV 
(see Fig. |7), constraining the cutoff energy for the EGRET 
power law with F = 2.38 to be E c < 14 GeV. Additionally, 
the H.E.S.S. limits significantly constrain the parameter space 
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Fig. 7. H.E.S.S. energy flux limits (99% c.l.) for pulsed emis- 
sion of the Vela pulsar (see Fig. [6] for point descriptions). 
The cross and open square deno te limits from the Du rham 
Mark 6 telescope (99.865% c.l., IChadwick etall l2000b. an d 
the CANGAROO detector (95% c.l.. lYoshikoshi et al.lll997b . 
respectively. The op en triangles resem ble the EGRET phase- 



averaged spectrum (IFierro et al.l I1998I) . The polar cap and 



oute r gap model cu r ves ar e generated according to the model 



of iRudak & Dvksl I l999. black dotted curve) and from 



Hiro tani et al.l (I2003L dashed curve), respectiv ely. The dotte d 



grey outer gap model curve is taken from lTakata et al.l (2006). 



for th e IC component of the outer gap model dHirotani et al 



2003), which was com puted with a sligh t change of the gap 
size D ± = 0.177 (see iHirotani et al.ll2003l for details) to bet- 
ter match_fheF^RETdata. Note that in the model calculation 
of Hirotani et al. (2003), no contribution from the photons of 
the optical puls e to the outer gap seed photons was considered 
(Hirotani 2006) which would result in a further increased IC 
component. On the other hand, a more detailed treatment of the 
outer gap aniso tropy would lead to a decreased IC contribution 
( Hirotanil2006h . A newer, two-dimensional version of the same 
outer gap model, which includes the contribution from the op- 
tical a nd infrared pulse, predicts a lower IC flux in the TeV 
range dTakata et al. 2006). The polar cap model again predicts 
pulsed emission only at energies far below the energy threshold 
of current ground-based y-ray detectors. 

PSRB1706-44. For this pulsar, limits on pulsed y-ra y 
emission were previously reported by IChadwick et al. (2000). 
Although H.E.S.S. data were mostly limited to 2-telescope data 
from the commissioning phase to which only the on-off-pulse 
method has been applied, more restrictive limits could be set 
(see Fig. [HJ, constraining the exponential cutoff energy to be 
E c < 7 1 GeV using the power law fit of the EGRET spectrum 
with T = 2.25. However, the IC component in the outer gap 
scenario of IHirotani et al.l d2003l) is left unconstrained by sev- 
eral orders of magnitude. 

PSRB1509-58. The rather strong unpulsed VHE y-ray 
emission of the PWN surrounding PSRB 1509-58 decreases 
the flux sensitivity for pulsed emission in comparison with 
other pulsars. The limits shown in Fig.|9]represent the first sen- 
sitive search for sub-TeV y-rays from this pulsar. The IC com- 
ponent of an outer gap model prediction (IHirotani & Shi bata 
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Fig. 8. H.E.S.S. energy flux limits (99% c.l.) for pulsed 
emission of PSRB 1706-44 for the on-off-pulse method (full 
points). The c ross denotes a flux Um it from the Durham group 
(99.865% c.l. JChadwick et alj|200fj|) . The op en triangles show 
the E GRET phase-averaged y-ray spectrum ( [Thompson et alj 
119961) . Polar cap a nd outer gap model c urves were calcu- 
lated according to IRudak & Dvksl 
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Fig. 9. H.E.S.S. energy flux limits (99% c.l.) for pulsed emis- 
sion of PSRB 1509-58 (see Fig.|6]for point descriptions). The 
open tria ngles resemble the EG RET upper limits for pulsed 
emission dThompson et al. 19991). The outer gap m odel curve 
(dotted) originates from Hirotani & Shibata] ( 200 lbl) . Fig. 6. 



2001b, dotted curve, parameters j\ - 0.6,^2 = 0) can be 
ruled out by the H.E.S.S. limits, although the curvature radi- 
ation component seems to be already marginally constrained 
by the EGRET limits at energies < 200 MeV. 

The remaining pulsars have not been detected in the MeV- 
GeV range and no detailed calculations based on any high- 
energy emission models are available for the TeV range. 

Independent of any y-ray emission scenario, for all selected 
pulsars the flux limits obtained imply an upper limit on the ratio 
of the total pulsed VHE y-ray luminosity compared to the pul- 
sars' spin-down power in the range 10~ 6 to 10~ 4 . This limit on 
the y-ray efficiency is significantly lower than the correspond- 
ing EGRET sensitivity in the GeV range. 
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Fig. 10. Differential upper limits on the energy flux of pulsed VHE y-ray emission for non-y-ray pulsars as measured by H.E.S.S. 
(at 99% confidence, point descriptions as in Fig. [6]). 



7. Conclusions 

H.E.S.S. data on young and energetic pulsars, including four of 
the eight known Ge V y-ray pulsars, have been searched for ev- 
idence of pulsed y-ray emission at energies between 100 GeV 
and lOOTeV. No conclusive evidence for pulsed emission has 
been found and differential upper limits on the pulsed flux were 
derived using a variety of different methods, constraining the 
pulsed flux for a wide range of possible pulse shapes and spec- 
tra in the VHE y-ray range. 

In the case of the well-studied y-ray pulsars Crab, Vela, 
and PSR B1509-58, the pulsed flux upper limits reported here 
significantly constrain the inverse Compton component of se- 
lected outer gap models for flux predictions in the TeV range 
which were available at the time of writing. However, the re- 
sults cannot rule out the outer gap scenario in general, since not 
all of the different calculations published give predictions for 
the IC component. Also, free model parameters such as seed 
photon density, inclination, gap geometry or viewing angles 
can be adjusted for the model cases considered here. However, 
it remains an open challenge for theory to provide reasonable 
parameters consistent with the flux limits reported here. 

For the complete sample of pulsars, containing a large frac- 
tion of the most energetic pulsars currently known, the absence 
of pulsed VHE y-ray emission already disfavours a significant 
contribution of the IC component to the energy loss mechanism 



of these pulsars. This sample exhibits great variety in different 
pulsar characteristics, providing valuable input to further stud- 
ies of the outer gap IC component. For example, in the case of 
the binary pulsar PSRB1259-63, a large spin-down luminos- 
ity is accompanied by a rather weak magnetic field, i.e. reduced 
y-ray absorption by pair creation, and additionally the photon 
field originating from the companion star provides a high den- 
sity target with known energy spectrum for a possible pulsed 
IC component, fixing an otherwise uncertain model parameter. 

Although in several cases there is spatial coincidence with 
extended TeV y-ray emission, pulsed emission is not detected 
in VHE y-rays. In particular, the flux upper limits derived are 
of the order of 10~ 4 to 10~ 6 of the pulsar spin-down flux, un- 
derlining the non-magnetospheric origin of the TeV radiation 
component and supporting the widely accepted scenario of an 
effective energy transport mechanism to, and strong particle ac- 
celeration in, the pulsar wind nebula. 

Finally, the upper limits given here imply a steep turnover 
of the pulsed high-energy spectrum at energies of a few tens of 
GeV for all pulsars studied which have established high-energy 
(> keV) emission, confirming a prediction common to all pul- 
sar models. Moreover, for the remaining pulsars, a turnover 
at sub-TeV energies would be likely even if pulsed GeV y- 
ray emission with spectral properties similar to that of the 
establishe d y-ray pulsars is detected with future sensitive in- 
struments (Thompson 2004). However, the different scenarios 
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Fig. 11. Differential upper limits at 99% confidence on 
the energy flux of unpulsed VHE y-ray emission for 
PSRJ1747-2958, PSR J 180 1-2451 and PSR J181 1-1925 
(from top to bottom) as measured by H.E.S.S. 

for pulsar emission differ significantly in their predictions of 
the exact shape and energy of the turnover. Thus, the search 
for pulsed y-ray emission from pulsars provides interesting 
prospects for future satellite-based and low-threshold ground- 
based y-ray instruments. 
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